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1.  SUMMARY 

In  months  7-12,  detailed  force  measurements  on  a  plunging  rectangular  wing  with  semi  aspect 
ratio  of  2,  at  low  Reynolds  numbers  of  1x104-3x104  have  been  conducted  in  a  water  tunnel  at 
various  non-dimensional  frequencies  and  amplitudes,  alongside  particle  image  velocimetry  (PIV) 
measurements  of  the  flow  at  the  midspan  plane  and  hotfilm  measurements  in  the  wake.  The 
presence  of  multiple  peaks  in  lift  has  been  identified  for  this  three-dimensional  wing,  thought  to 
be  related  to  the  natural  shedding  frequency  of  the  stationary  wing.  Wing/vortex  and 
vortex/vortex  interactions  have  been  identified  which  may  also  contribute  to  the  selection  of 
optimal  frequencies.  Lift  enhancement  is  observed  to  become  more  notable  with  increasing 
amplitude,  coupled  with  a  slight  shift  in  the  peak  frequencies,  to  lower  Strouhal  numbers,  with 
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increasing  angle  of  attack.  Despite  the  highly  three-dimensional  nature  of  the  flow,  lift 
enhancements  up  to  1 80%  are  possible. 

The  unsteady  aerodynamics  of  low  aspect  ratio  membrane  wings  (rectangular  wing  and  delta 
wing)  has  been  investigated  in  a  Reynolds  number  range  of  2.4*  104  -  5.9X104.  Membrane 
deformations  were  measured  by  using  a  Digital  Image  Correlation  (DIC)  system.  Force 
measurements  were  conducted  for  membrane  wings  and  compared  with  those  of  rigid  wings  of 
similar  geometries.  The  comparison  indicates  that  the  effect  of  flexibility  benefits  the  rectangular 
wing  more  than  the  delta  wing  by  increasing  the  maximum  normal  force  as  well  as  the  slope  by  a 
larger  amount.  Oscillations  of  the  membrane,  mode  shapes,  and  the  corresponding  vibration 
frequencies  were  investigated.  The  combination  of  tip  vortices  and  leading-edge  vortex  shedding 
results  in  a  mixture  of  streamwise  and  spanwise  vibrational  modes.  The  Strouhal  numbers  of 
membrane  vibration  frequencies  are  on  the  order  of  unity,  and  have  similar  range  to  those  of 
two-dimensional  airfoils.  The  results  suggest  possible  coupling  with  wake  instabilities.  The 
Particle  Image  Velocimetry  (PIV)  measurements  indicate  stronger  vortices  for  membrane  wings, 
which  contribute  significantly  to  the  total  lift.  The  experimental  data  obtained  in  the  past  6 
months  are  being  used  for  comparisons  with  the  computational  simulations  carried  out  by  Dr. 
Visbal’s  group. 

2.  EXPERIMENTAL  APPARATUS  AND  METHODS 

2.1  Plunging  Rigid  Low  Aspect  Ratio  Wings 

Experiments  were  carried  out  in  a  free-surface  closed-loop  water  tunnel  (Eidetics®  Model 
1520)  at  the  University  of  Bath.  The  forced  plunging  motion  of  a  rectangular  wing  (semi  aspect 
ratio  sAR=2)  with  NACA0012  cross-section  was  provided  by  a  ‘shaker  mechanism’,  consisting 
of  an  AC  0.3 7kW  Motovario  three-phase  motor,  using  a  5:1  gearbox  reduction  and  a  basic  rotary 
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to  linear  mechanism.  An  illustration  of  the  rig  is  provided  in  Figure  1.  Flow  measurements  were 
obtained  using  a  Digital  Particle  Image  Velocimetry  (DPIV)  system.  The  TSI  Inc.  2D-PIV 
system  comprises  of  a  dual  ND:YAG  120  mJ  pulsed  laser,  and  was  used  to  illuminate  a 
streamwise  plane.  The  flow  was  seeded  with  8-12  pm  hollow  glass  particles  and  the  illuminated 
plane  was  then  captured  by  a  2MP  Powerview  Plus  12bit  CCD  camera.  A  rotary  encoder  was 
used  to  send  a  pulse  360  times  a  second  each  with  a  constant  time  interval.  The  combination  of 
an  trigger  system,  to  isolate  a  desired  pulse  from  the  rotary  encoder  and  the  LaserPulse 
synchroniser,  was  used  to  fire  the  laser  at  the  required  point  in  the  cycle.  Both  phase-averaged 
and  time-averaged  measurements  were  conducted.  Captured  images  were  then  processed  using  a 
commercially  available  software,  TSI  Insight  3G,  to  produce  a  series  of  velocity  vectors.  Force 
measurements  were  conducted  using  a  two  component  binocular  load  cell.  An  AD/DA  converter 
converts  the  signals  before  being  sent  to  the  data  acquisition  card. 

2.2  Low  Aspect  Ratio  Membrane  Wings 

The  experiments  were  carried  out  in  a  low-speed,  closed-loop  open-jet  wind  tunnel  with  a 
circular  working  section  of  760  mm  in  diameter,  located  at  the  Department  of  Mechanical 
Engineering  of  University  of  Bath.  After  discussions  with  AFRL  collaborators  (Miguel  Visbal), 
two  wings  were  tested,  a  rectangular  wing  with  aspect  ratio  of  AR  =  2  (a  chord  length  of  c  =  68.8 
mm)  and  a  delta  wing  with  sweep  angle  of  A  =  50°  (a  chord  length  of  c  =  82.8  mm).  The 
dimensions  are  shown  in  Figures  2a  and  2b.  Each  wing  was  composed  of  a  rigid  frame  and  an 
attached  membrane.  The  wing  was  attached  to  a  shaft  through  a  strain  gauge  force  balance  for 
making  direct  force  measurements  in  the  normal  directions.  For  the  purpose  of  comparison,  force 
measurements  were  also  conducted  on  rigid  wings.  The  shaft  was  driven  by  a  computer 
controlled  step-motor  for  changing  the  wing  incidences.  A  high  frame  rate  Digital  Image 
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Correlation  (DIC)  system,  which  employs  two  pre-calibrated  CCD  cameras,  was  used  to  measure 
the  time-accurate  membrane  displacements.  The  DIC  system  employs  a  mathematical  correlation 
method  to  measure  deformation  on  an  object's  surface  [1-3].  Experiments  were  conducted  over  a 
range  of  angle  of  attack  and  various  Reynolds  numbers.  Quantitative  flow  measurements  were 
undertaken  using  Digital  Particle  Image  Velocimetry  (DPIV).  The  laser  sheet  was  placed  parallel 
to  the  freestream  flow  at  the  mid-span  location  for  the  measurements  in  a  streamwise  plane  as 
shown  in  Figure  3a,  and  perpendicular  to  the  freestream  very  close  to  the  trailing  edge  (at  2  mm 
distance)  for  the  cross-flow  measurements  as  shown  in  Figure  3b.  The  commercial  software  TSI 
Insight  v6.0  and  a  fast  Fourier  transform  (FFT)  cross-correlation  algorithm  were  used  to  analyze 
the  captured  images. 

3.  MAIN  FINDINGS 

3.1  Plunging  Rigid  Low  Aspect  Ratio  Rectangular  Rigid  Wing 

Figure  4  presents  the  time-averaged  lift  and  drag  of  the  wing  oscillating  at  an  amplitude  of 
h=0.025  and  an  incidence  of  a=20°  and  indicated  Reynolds  numbers.  A  substantial  increase  in 
lift  is  recorded  with  dominant  peaks  occurring  at  Stc~0.45-0.5  and  Stc~0.9.  The  corresponding 
lift  enhancements  are  67%  and  75%,  respectively,  over  the  stationary  case,  at  Re=lx104. 
Increasing  the  frequency  beyond  the  second  peak,  results  in  a  plateau  being  reached,  albeit  with 
underlying  attenuated  undulations.  Within  a  range  of  0<Stc<1.2  the  effect  of  varying  the 
operating  Reynolds  number  between  Re=lx104  and  Re=3x104  has  little  effect,  suggesting  that 
the  phenomenon  is  aerodynamic  and  not  related  to  any  structural  resonance  or  setup-related 
issues.  Although  not  primarily  the  motivation  for  the  current  study,  drag  measurements  have  also 
been  presented,  given  that  there  exists  considerable  interest  in  the  propulsion  of  oscillating 
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wings.  A  peak  in  drag  occurs  simultaneously  with  the  initial  peak  in  lift,  followed  by  a  gradual 
reduction  with  increased  frequency,  even  so,  no  thrust  is  produced  within  the  frequency  range 
tested. 

The  phase-averaged  vorticity  and  velocity  magnitude,  at  the  midspan  plane,  are  presented  in 
Figure  5a  and  Figure  5b,  respectively.  Figure  5  illustrates  the  presence  of  multiple  vortices 
appearing  on  the  suction  surface  of  the  wing,  the  number  of  which  depends  on  the  frequency  of 
the  oscillation.  Multiple  distinct  LEVs  have  been  observed  in  two-dimensional  flows  [4,5],  and 
evidently  still  present  to  at  least  one  chord  length  outboard  of  the  root  of  this  low  aspect  ratio 
wing.  At  the  location  of  the  initial  peak  in  lift,  Stc=0.45,  one  can  observe  the  size  and  location  of 
the  generated  vortex  (LEVi).  Interestingly,  at  this  frequency,  a  single  vortex  is  observed  over  the 
wing,  in  comparison  to  the  two  visible  leading-edge  vortices  over  the  wing  at  the  second  peak  in 
lift  (Stc=0.9).  At  subsequent  multiples  of  Stc=0.9,  an  additional  leading  edge  vortex  appears, 
given  that  the  increased  forcing  frequency  reduces,  by  a  similar  factor,  the  distance  travelled  per 
cycle.  This  seems  to  suggest  that  these  optimal  frequencies  are  selected  based  on  the  number  of 
vortices  captured  by  the  wing.  A  coupling  between  the  wake  instabilities  and  vortex/wing 
interactions  may  therefore  determine  these  optimal  frequencies  for  lift. 

3.1.1  Effect  of  oscillation  frequency 

A  previous  study  [6]  on  the  aerodynamic  performance  of  an  oscillating  two-dimensional 
NACA0012  profile,  elicited  the  existence  of  several  frequencies  at  which  the  time-averaged  lift 
is  locally  maximised.  The  forcing  frequencies  corresponding  to  such  peaks  were  suggested  as 
being  sub-harmonic,  fundamental  and  harmonics  of  the  natural  vortex  shedding  frequency  in  the 
wake  of  the  stationary  wing  at  a=15°.  Figure  6  presents  a  comparison  of  the  time-averaged  lift 
performance  of  the  two-dimensional  airfoil,  to  the  current  finite  wing  of  sAR=2,  in  which  both 
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wings  are  subject  to  an  oscillating  amplitude  of  h=0.15,  set  to  a=20°,  the  latter  representative  of 
a  post-stall  condition.  Evidently,  optimal  frequencies  exist  for  the  finite  wing,  roughly  adhering 
to  those  of  the  airfoil.  A  drop  of  38-40%,  in  relation  to  the  local  maxima  in  lift  is  encountered. 
As  the  aspect  ratio  is  reduced,  an  increasingly  noteworthy  interaction  between  wing  tip  vortices 
and  leading  edge  vortices  occurs  [7].  If  we  consider  these  imparted  peaks  to  be  most  effective  in 
two-dimensional  flows,  added  unsteadiness  from  the  aforementioned  interaction  and  thus 
reduction  of  assimilative  two-dimensional  flow  neighbouring  the  midspan,  with  reduced  aspect 
ratio,  will  diminish  the  effectiveness  of  the  generated  vortices,  resulting  in  a  less  distinguishable 
peak.  This  is  evident  from  the  comparison  made  in  Figure  6,  through  the  apparent  attenuation 
and  graduality  of  the  peaks  of  the  low  aspect  ratio  wing. 

It  was  suggested  that  a  natural  vortex  shedding  frequency  of  Stc=0.83,  for  the  airfoil  case  at 
a=15°  [6],  and  its  sub-harmonic  and  first-harmonic  correspond  to  the  optimal  frequencies  in  lift. 
Hot-film  measurements  were  performed  to  investigate  wake  instabilities  for  the  current  wing.  A 
typical  signal  for  the  presented  case  of  a=20°  and  Re=lxl04  is  shown  in  Figure  7.  The  quasi- 
periodic  nature  of  the  flow  at  the  midspan,  is  akin  to  the  supercritical  mode  observed  by  Huang 
and  Fin  [8]  and  adequate  enough  to  produce  a  distinct  spectral  peak,  presented  in  Figure  8.  The 
measured  frequencies  at  three  different  angles  attack  are  presented  in  Figure  9  as  hollow 
triangular  symbols.  A  comparison  has  also  been  made  within  Figure  9  with  the  vortex  shedding 
frequencies  reported  by  Huang  and  Fin  [8]  using  the  same  NACA0012  profile  at  Re=l.lxl04, 
albeit  with  a  larger  semi  aspect  ratio  of  5.  The  results  reveal  a  surprisingly  remarkable  fit,  despite 
the  lower  aspect  ratio  of  the  current  wing.  Superimposed  within  the  same  graph  are  the  locations 
of  the  first  and  second  time-averaged  lift  peaks  for  the  three  angles  of  attack  that  have  been 
tested  in  the  present  study.  Figure  9  shows  that  the  first  peak  for  a=20°,  is  very  close  to  the 
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natural  vortex  shedding  frequency  in  the  wake.  However,  for  a=10°,  the  data  seem  to  suggest 
that  the  first  peak  could  be  the  sub-harmonic  of  the  natural  vortex  shedding  frequency.  The 
interpretation  is  even  more  difficult  for  a=15°,  for  which  the  first  peak  corresponds  to  the  sub 
harmonic  in  the  two-dimensional  airfoil  case  [6],  It  is  apparent  that  the  peaks  in  lift  are  in  the 
same  order  as  the  frequencies  of  the  wake  instabilities  of  the  stationary  wing.  However,  the  first 
peak  remains  remarkably  unchanged  for  all  three  incidences,  while  the  fundamental  frequency  of 
the  wake  instability  varies. 

3.1.2  Effect  of  amplitude 

Figure  10  presents  lift  and  drag  measurements  for  the  wing  oscillating  with  an  amplitude  of 
h=0.15  for  all  Re  cases.  This  is  presented  with  the  measurements  for  h=0.025  operating  at 
Re=lx104.  There  is  an  immediate  benefit  to  operating  at  this  higher  amplitude.  Quantitatively, 
what  was  originally  observed  to  be  a  67%  increase  in  lift  at  Stc=0.45  for  h=0.025,  transforms  to 
an  increase  of  126%  over  the  static  case  for  h=0.15.  Moreover,  the  change  in  amplitude  has, 
within  experimental  uncertainty,  unaffected  the  required  forcing  frequency  of  the  local  peak  in 
lift.  Again,  little  difference  can  be  inferred  between  the  Reynolds  numbers  tested.  The  switch 
from  drag  to  thrust  occurs  within  the  operational  range  of  Strouhal  numbers,  at  Stc~l .  I  - 1 .25. 

Figure  11  shows  the  phase-averaged  vorticity  contour  plots  for  varying  Stc,  at  various 
locations  in  the  cycle.  Analogously  to  the  lower  amplitude  case,  as  the  frequency  is  increased, 
the  vortices  form  closer  to  the  surface  with  greater  coherency.  Notice  the  LEV  formed  at  Stc=0.2 
compared  to  Stc=0.45.  The  increased  frequency  evokes  a  relatively  lower  convection  distance 
within  the  cycle,  allowing  the  vortex  to  remain  close  to  the  upper  surface,  through  to  the  bottom 
of  the  cycle.  The  delayed  shedding  of  the  vortex  would  imply  a  more  effective  form  of  lift 
production,  resulting  in  the  monotonic  rise  in  lift  up  to  a  Strouhal  number  of  Stc=0.45.  Further 


increasing  the  forcing  frequency,  encourages  this  delay,  eventually  preventing  the  leading  edge 
vortex  from  passing  the  trailing  edge  vortex  formed  in  the  following  cycle.  Once  again,  in  a 
similar  fashion  to  h=0.025,  two  LEVs  are  present  on  the  upper  surface,  at  the  bottom  of  the 
oscillation,  for  Stc=0.9  (Figure  11a).  Although  the  LEV,  from  the  preceding  cycle,  is  located 
further  away  from  the  upper  surface  of  the  wing,  than  that  of  h=0.025,  the  streamwise  convection 
distance  is  very  similar.  A  further  increase  in  frequency  has  the  effect  of  causing  a  destructive 
interaction  between  the  wing  and  LEV,  similar  to  that  of  the  two-dimensional  airfoil.  This  is 
clearly  observed  for  Stc=2.0,  for  which  the  LEV  is  absent  at  the  top  of  the  oscillation  (t/T=1.0). 
This  in  effect  precludes  any  further  formation  of  multiple  vortices  on  the  upper  surface. 

3.2  Unsteady  Aerodynamics  of  Low  Aspect  Ratio  Membrane  Wings 

Figure  12  presents  the  variation  of  the  maximum  magnitudes  of  z-displacement  as  a  function 
of  angle  of  attack  of  the  rectangular  wing  for  different  freestream  velocities.  For  the  lowest 
freestream  velocity  Uoo  =  5  m/s,  the  maximum  displacement  increases  gradually  up  to  about  a  = 
20°,  then  stays  roughly  constant.  The  membrane  generates  a  larger  magnitude  of  the  maximum 
camber  as  the  freestream  velocity  is  increased.  Notably,  there  is  also  camber  at  zero  angle  of 
attack  particularly  at  U„  =10  m/s.  The  variation  of  the  maximum  magnitudes  of  the  normalized 
displacement  of  the  delta  wing  as  a  function  of  angle  of  attack  is  shown  in  Figure  13  for  different 
freestream  velocities.  Generally,  the  delta  wing  exhibits  similar  trend  of  the  variations  in 
displacements  to  those  of  the  rectangular  wing,  but  with  smaller  magnitude.  For  U,x  =  5  m/s  and 
at  small  incidences,  the  slope  is  relatively  steep  which  means  that  the  membrane  cambering  is 
more  sensitive  to  the  change  in  angles  of  attack  within  this  range.  For  this  freestream  velocity, 
the  maximum  displacement  increases  gradually  up  to  about  a  =  21°.  Once  again  the  membrane 
shows  a  larger  magnitude  of  the  maximum  displacement  as  the  freestream  velocity  is  increased. 
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However,  the  rate  of  the  increase  of  the  maximum  displacement  zmax/c  with  angle  of  attack  a  is 
smaller  for  the  higher  freestream  velocities. 

The  variation  of  the  time-averaged  normal  force  coefficient  is  shown  as  a  function  of  angle  of 
attack  in  Figure  14  for  the  membrane  and  rigid  rectangular  wings  at  different  freestream 
velocities.  Even  though  the  displacement  slopes  (Figure  12)  are  not  linear,  the  normal  force  slope 
is  nearly  linear  until  the  stall.  The  flexible  membrane  wings  have  higher  slopes  than  the  rigid 
wing,  and  the  stall  angles  are  only  slightly  affected.  These  favorable  larger  slopes  are  not 
unexpected  as  the  induced  wing  camber  increases  with  increasing  angle  of  attack.  However,  the 
corresponding  increase  in  the  camber  alone  does  not  explain  the  significant  increase  in  the  force 
slope.  In  addition,  it  is  observed  that  the  membrane  wings  start  to  generate  the  lift  at  a  =  0°.  This 
is  consistent  with  the  wing’s  maximum  displacement  (Figure  12)  which  shows  a  small  camber  at 
zero  angle  of  attack.  At  the  highest  freestream  velocity,  Uoo  =10  m/s,  the  membrane  wing 
reaches  the  highest  maximum  normal  force  coefficient  Cn>max  of  approximately  1.26,  which  is 
about  1 .5  times  of  that  of  the  rigid  wing.  Similar  to  the  rectangular  wing,  the  flexible  delta  wing 
exhibits  a  larger  time-averaged  normal  force  coefficient  than  the  rigid  delta  wing  as  shown  in 
Figure  15.  At  the  highest  freestream  velocity,  Uoo  =10  m/s,  the  membrane  wing  reaches  the 
highest  Cn  max  of  1.16,  which  is  approximately  1.2  times  that  of  the  rigid  wing. 

Figure  16  shows  the  frequency  spectra  computed  from  the  z-displacement  of  the  membrane 
sampled  at  the  point  of  maximum  standard  deviation  at  Uoo  =10  m/s.  For  small  incidences,  the 
dominant  normalized  frequencies  (Strouhal  numbers  Sr  =  fc/Uoo)  are  rather  difficult  to  identify, 
suggesting  that  the  magnitudes  of  vibration  are  very  small.  However,  starting  around  a  =  7°  -  9°, 
the  dominant  normalized  frequencies  fall  in  a  well-defined  band,  and  on  the  order  of  unity. 
Within  certain  ranges  of  angles  of  attack,  there  seems  to  be  higher  modes.  Of  course,  the  change 
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in  the  dominant  frequencies  is  related  to  the  change  in  the  mode  shapes.  At  high  incidences, 
where  the  second  mode  is  always  observed  for  all  the  velocities  tested,  the  Strouhal  numbers  are 
in  a  similar  range  with  those  of  the  two-dimensional  membrane  airfoils  [9].  For  two-dimensional 
membrane  airfoils,  the  oscillation  in  this  range  suggests  the  coupling  with  the  vortex  shedding  in 
the  wake.  Therefore,  the  present  results  suggest  a  possible  coupling  of  membrane  oscillation  for 
the  rectangular  wing  with  the  wake  instabilities. 

Figure  17  shows  the  comparison  of  the  magnitude  of  the  time-averaged  velocity  field 
superimposed  on  the  streamline  pattern  for  the  rectangular  wings  at  selected  angles  of  attack. 
The  data  were  taken  for  the  rigid  wing  at  Uoo  =  5  m/s,  and  the  membrane  wing  at  Uoo  =  5  m/s  and 
10  m/s.  Starting  with  small  incidence  a  =  5°,  the  flow  remains  attached  for  all  wings.  At  higher 
angle  of  attack  a  =  10°,  the  flow  on  the  rigid  wing  becomes  separated.  For  a  =  16°,  the  shear 
layer  moves  away  from  the  surface  for  the  rigid  wing,  while  the  flow  on  the  membrane  wing  is 
separated  but  the  shear  layer  remains  close  to  the  surface.  It  is  also  seen  that  at  this  angle  of 
attack,  with  increasing  freestream  velocity  to  U,x  =  10  m/s,  the  shear  layer  for  the  membrane 
wing  is  getting  closer  to  the  wing  surface  while  the  camber  of  the  wing  increases.  The  effect  of 
increasing  freestream  velocity  is  consistent  with  the  previous  work  for  two-dimensional 
membrane  airfoils  [9,10].  The  separated  shear  layer  moves  further  away  from  the  surface  with 
increasing  angle  of  attack  to  a  =  23°  for  all  wings.  Again,  the  separation  region  is  largest  for  the 
rigid  wing,  and  smallest  for  the  membrane  wing  with  Uoo  =10  m/s.  For  the  latter,  the  closed 
streamline  pattern  (recirculation  region)  is  seen  over  the  wing  surface,  while  it  is  located  further 
downstream  for  the  lower  freestream  velocity,  and  even  further  away  from  the  trailing  edge  for 
the  rigid  wing. 
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Figure  18  shows  the  magnitude  of  the  time-averaged  velocity  in  a  cross-flow  plane  near  the 
trailing  edge  for  the  rectangular  rigid  and  membrane  wings.  The  cross-flow  velocity  is  seen  to 
increase  with  increasing  flexibility.  This  suggests  stronger  tip  vortices  with  increasing  camber. 
The  variation  of  normalized  circulation  of  tip  vortices  (not  shown),  calculated  from  the  area 
integral  of  vorticity  distributions  in  the  cross-flow  plane  based  on  the  PIV  measurements, 
suggested  that  membrane  wing  has  a  larger  magnitude  of  the  circulation  than  the  rigid  wing  for 
both  freestream  velocities.  Since  the  tip  vortices  (and  vortex  lift)  contribute  significantly  to  the 
total  lift,  this  may  be  an  additional  mechanism  for  the  increased  lift  with  increasing  membrane 
deformation.  Whereas  the  normalized  circulation  of  the  rigid  wing  remains  virtually  unchanged 
with  increasing  freestream  velocity  (Reynolds  number),  the  normalized  circulation  increases  for 
the  membrane  wing. 
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Figure  1.  An  illustration  of  the  water  tunnel  experimental  setup  and  apparatus 


Figure  2.  Schematic  of  the  wings  with  cross-section  of  the  rigid  frame;  a)  Rectangular  wing;  b) 
Delta  wing. 


a)  b) 

Figure  3. Schematic  of  PIY  measurements  for:  a)  streamwise  plane  at  a  mid-span  and  b) 
crossflow  plane  near  trailing-edge. 
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Figure  4.  Time  averaged  lift  (left)  and  drag  (right)  measurements  for  the  wing  oscillating  at 
h=0.025,  a=20°  and  indicated  Re. 


a)  b) 


Figure  5.  Phase-averaged  a)  vorticity  contours  and  b)  streamlines  for  h=0.025,  Re=lx104  and 
a=20°  at  t/T=0.5  (bottom  of  the  motion).  Each  LEVi  represents  a  vortex  formed  in  a  separate 
cycle. 
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Figure  6.  A  comparison  between  the  time- 
averaged  lift  coefficients  of  the  2D  and  sAR=2 
NACA0012  wing  oscillating  at  h=0.15  and 
a=20°,  for  Re=lxl04. 
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Figure  8.  Spectral  analysis  of  the  wake 
velocity  for  the  stationary  wing  at  a=20°. 
x/c=0.5,  y/c=-0.2  and  z/c=l.  The 

streamwise  (x)  and  cross-stream  (y)  distance 
are  measured  from  the  trailing  edge. 
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Figure  7.  Time-series  of  the  wake  velocity  for 
the  stationary  wing  at  a=20°.  x/c=0.5,  y/c=- 
0.2  and  z/c=l.  The  streamwise  (x)  and  cross¬ 
stream  (y)  distance  are  measured  from  the 
trailing  edge. 
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Figure  9.  Variation  of  peaks  in  lift  coefficient 
for  h=0.15,  Re=l><104,  in  relation  to  the 
natural  vortex  shedding  frequencies  of  the 
stationary  wing. 
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Figure  10.  Time-averaged  lift  (left)  and  drag  (right)  measurements  for  h=0.15,  a=20°  and 
indicated  Re. 
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Figure  11.  Phase-averaged  normalised  vorticity  plots  for  h=0.15,  a=20°,  Re=lx104  and  stated 
Stc,  at  a)  t/T=0.25,  b)  t/T=0.5,  c)  t/T=0.75,  d)  t/T=1.0. 
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Figure  12.  Maximum  magnitude  of  time-  Figure  13.  Maximum  magnitude  of  time- 
averaged  displacement  as  a  function  of  averaged  displacement  as  a  function  of 
incidence  for  the  rectangular  wing.  incidence  for  the  delta  wing. 


Figure  14.  Normal  force  coefficient  as  a  Figure  15.  Normal  force  coefficient  as  a 
function  of  incidence  for  the  flexible  and  rigid  function  of  incidence  for  the  flexible  and  rigid 
rectangular  wings.  delta  wings. 


Figure  16.  Power  spectra  of  the  membrane  oscillations  as  a  function  of  angle  of  attack  for  the 
rectangular  wing  at  Uoo  =10  m/s. 
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a) 


b) 


c) 


Figure  17.  Magnitude  of  the  time-averaged  velocity  and  streamlines  in  a  mid-span  plane  for  the 
rectangular  wings;  a)  rigid  wing  U.n  =  5  m/s;  b)  membrane  wing  UB  =  5  m/s;  c)  membrane  wing 
Uoo  =  10  m/s.  Flow  is  from  right  to  left. 
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Figure  18.  Magnitude  of  the  time-averaged  velocity  in  a  crossflow  plane  near  the  trailing  edge 
for  the  rectangular  wings;  a)  rigid  wing  Uoo  =  5  m/s;  b)  membrane  wing  U«,  =  5  m/s;  c)  membrane 
wing  Uoo  =  10  m/s. 
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